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Abstract 



' Quantum adiabatic processes -that keep constant the populations in the instan- 

taneous eigenbasis of a time-dependent Hamiltonian- are very useful to prepare and 
manipulate states, but take typically a long time. This is often problematic because 
^ ■ decoherence and noise may spoil the desired final state, or because some applications 

. require many repetitions. "Shortcuts to adiabaticity" are alternative fast processes 



which reproduce the same final populations, or even the same final state, as the 
adiabatic process in a finite, shorter time. Since adiabatic processes are ubiquitous, 
the shortcuts span a broad range of applications in atomic, molecular and optical 
physics, such as fast transport of ions or neutral atoms, internal population control 
and state preparation (for nuclear magnetic resonance or quantum information), 
cold atom expansions and other manipulations, cooling cycles, wavepacket split- 
ting, and many-body state engineering or correlations microscopy. Shortcuts are 
also relevant to clarify fundamental questions such as a precise quantification of the 
third principle of thermodynamics and quantum speed limits. We review different 
theoretical techniques proposed to engineer the shortcuts, the experimental results, 
and the prospects. 
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quantum state engineering, transport engineering of cold atoms, ions, and 
Bose-Einstein condensates, wave packet splitting, third principle of 
thermodynamics, transitionless tracking algorithm, fast expansions 
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1 Introduction 

The expression "sliortcuts to adiabaticity" (STA) was recently introduced in 
Chen et al. (2010b), to describe protocols that speed up a quantum adiabatic 
process, usually, although not necessarily, through a non-adiabatic route. ^ 
There the Lewis-Riesenfeld invariants were used to inverse engineer the time 
dependence of a harmonic oscillator frequency between predetermined initial 
and final values so as to avoid final excitations. That paper and its companion 
on Bose- Einstein condensates (Muga et al., 2009) have indeed triggered a surge 
of activity, not only for harmonic expansions (Chen and Muga, 2010; Muga 
et al., 2010; Stefanatos et al, 2010; Schaff et al., 2010, 2011a; del Campo, 
2011a; Schaff et al., 2011b; Stefanatos et al., 2011; Torrontegui et al., 2012c; 
Fasihi et al., 2012; Torrontegui et al., 2012a; del Campo and Boshier, 2012; Ste- 
fanatos and Li, 2012), but for atom transport (Torrontegui et al., 2011, 2012d; 
Chen et al., 2011b; Bowler et al., 2012), quantum computing (Sarandy et al., 
2011), quantum simulations (Lau and James, 2012), optical lattice expansions 
(Yuce, 2012; Ozcakmakli and Yuce, 2012), wavepacket splitting (Torrontegui 
et al., 2012b), internal state control (Chen et al., 2011a; Ibafiez et al., 2011; 
Ruschhaupt et al., 2012; Ban et al., 2012; Ibafiez et al., 2012a; Giingordii et 
al., 2012), many-body state engineering (del Campo, 2011b; del Campo and 
Boshier, 2012; del Campo et al., 2012; Julia-Diaz et al., 2012), and other ap- 
plications such as sympathetic cooling of atomic mixtures (Choi, Onofrio and 
Sundaram, 2011; Choi et al., 2012), or cooling of nanomechanical resonators 
(Li et al., 2011; Zhang et al., 2012a). In fact several works had previously 
or simultaneously considered to speed up adiabatic processes making use of 
different techniques. For example Demirplak and Rice (2003, 2005, 2008) and 
Berry (2009) proposed the addition of counterdiabatic terms to a reference 
Hamiltonian Ho to achieve adiabatic dynamics with respect to Hq. This "tran- 
sitionless tracking algorithm" (Berry, 2009) has been applied to manipulate 

^ The word "adiabatic" may have two different meanings: the thermodynamical one 
(no heat transfer between system and environment) and the quantum one, as stated 
by Born and Fock (1928) in the adiabatic theorem: "a physical system remains in its 
instantaneous eigenstate when a given perturbation is acting on it slowly enough and 
if there is a gap between the eigenvalue and the rest of the Hamiltonian's spectrum" . 
Here we shall always understand "adiabatic" in the quantum mechanical sense. 
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the populations of two-level systems (Demirplak and Rice, 2005; Berry, 2009; 
Chen et al., 2010c; Bason et al, 2012; Zhang et al., 2012b). Another tech- 
nique to design laser pulses for fast population transfer is parallel adiabatic 
passage (Gulerin et al., 2002; Vasilev et al., 2009; Dridi et al., 2009; Gulerin 
et al., 2011). Convert et al. (2008a) designed trap motions in order to perform 
non-adiabatic fast transport of atomic cold clouds. Also, Masuda and Naka- 
mura (2010) developed a "fast-forward technique" for several manipulations 
on wavepackets such as expansions, transport or splitting of Bose-Einstein 
condensates. Related work had also been carried out for wave packet splitting 
making use of optimal control (Hohenester et al., 2007; Grond et al., 2009a,b), 
and in the context of quantum refrigerators, to find fast "frictionless" expan- 
sions (Salamon et al., 2009; Rezek et al., 2009). For recent developments on 
this line stimulated by invariant-based engineering results see Kosloff and Feld- 
mann (2010); Hoffmann et al. (2011); Levy and Kosloff (2012); Feldmann and 
Kosloff (2012). 

The considerable number of publications on the subject, and a recent Confer- 
ence on "Shortcuts to adiabaticity" held in Bilbao (16-20 July 2012) demon- 
strate much current interest, not only within the cold atoms and atomic 
physics communities but also from fields such as semiconductor physics and 
spintronics (Ban et al., 2012). Indeed adiabatic processes are ubiquitous, so we 
may expect a broad range of applications, even beyond the quantum domain, 
since some of the concepts are easy to translate into optics (Lin et al., 2012; 
Tseng and Chen, 2012) or mechanics (Ibanez et al., 2011). Apart from the 
practical applications, the fundamental implications of shortcuts on quantum 
speed limits (Bender et al., 2007; Bason et al., 2012; Uzdin et al., 2012), time- 
energy uncertainty relations (Chen and Muga, 2010), multiple Schrodinger 
pictures (Ibanez et al., 2012a), and the quantification of the third principle of 
thermodynamics and of maximal cooling rates (Salamon et al., 2009; Rezek et 
al., 2009; Chen and Muga, 2010; Kosloff and Feldmann, 2010; Hoffmann et al., 
2011; Levy and Kosloff, 2012; Feldmann and Kosloff, 2012) are also intriguing 
and provide further motivation. 

In this review we shall first describe different approaches to STA in Sec. 2. 
While the main goal there is to construct new protocols for a fast manipula- 
tion of quantum states avoiding final excitations, additional conditions may 
be imposed. For example, ideally these protocols should not be state specific 
but work for an arbitrary state ^ . They should also be stable against per- 
turbations, and keep the values of the transient energy and other variables 
manageable throughout the whole process. Several applications are discussed 



^ Contrast this to the quantum brachistochrone (Bender et al., 2007), in which the 
aim is to find a time-independent Hamiltonian that takes a given initial state to 
a given final state in minimal time. Studies of "quantum speed limits" adopt in 
general this state-to-state approach, as in Uzdin et al. (2012). 
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in Sees. 3 to 6. We have kept a notation consistency within each Section but 
not throughout the whole review, following when possible notations close to 
the original publications. 



2 General Formalisms 



2.1 Invariant Based Inverse Engineering 



Lewis- Riesenf eld invariants. - The Lewis and Riesenfeld (1969) theory is ap- 
plicable to a quantum system that evolves with a time-dependent Hermitian 
Hamiltonian H{t), which supports a Hermitian dynamical invariant I{t) sat- 
isfying 

zh^-[H{t),I{t)]=0. (1) 

Therefore its expectation values for an arbitrary solution of the time- dependent 
Schrodinger equation ih-^\'^{t)) = H{t)\'^{t)), do not depend on time. I{t) 
can be used to expand |^E'(t)) as a superposition of "dynamical modes" |?/'„(t)), 

I^W) = ECnl^n(t)), \Mt)) = e^""W|0^(t)), (2) 
n 

where n = 0,1,...; Cn are time- independent amplitudes, and \(pn(t)) are or- 
thonormal eigenvectors of the invariant I{t), 

i{t) = Y.\Mt))K{Mt)l (3) 

n 

The A„ are real constants, and the Lewis- Riesenfeld phases are defined as 
(Lewis and Riesenfeld, 1969) 

«n(t) = ^ (^Mt')\^^^ - H{t')\Mt'))dt'. (4) 

We use for simplicity a notation for a discrete spectrum of I{t) but the general- 
ization to a continuum or mixed spectrum is straightforward. We also assume 
a non-degenerate spectrum. Non- Hermitian invariants and Hamiltonians have 
been considered for example in Gao et al. (1991, 1992); Lohe (009); Ibanez et 
al. (2011). 

Inverse engineering.- Suppose that we want to drive the system from an initial 
Hamiltonian H{0) to a final one H{tf), in such a way that the populations 
in the initial and final instantaneous bases are the same, but admitting tran- 
sitions at intermediate times. To inverse engineer a time-dependent Hamil- 
tonian H{t) and achieve this goal, we may first define the invariant through 
its eigenvalues and eigenvectors. The Lewis- Riesenfeld phases an{t) may also 
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be chosen as arbitrary functions to write down the time- dependent unitary 
evolution operator U, 

t/ = 5]e-"W|0.(t))(0„(O)|. (5) 

n 

U obeys ihU = H{t)U, where the dot means time derivative. Solving formally 
this equation for H{t) = ihUU\ we get 

H{t) = -hJ2 |0n(t))«n(0n(t)| + ^/^E l^t^n (0n(t) | • (6) 
n n 

According to Eq. (6), for a given invariant there are many possible Hamilto- 
nians corresponding to different choices of phase functions an{t). In general 
/(O) does not commute with H{0), so the eigenstates of /(O), |0„(O)), do not 
coincide with the eigenstates of H{0). H(tf) does not necessarily commute 
with I{tf) either. If we impose [/(0),i7(0)] = and [I{tf),H{tf)] = 0, the 
eigenstates will coincide, which guarantees a state transfer without final ex- 
citations. In typical applications the Hamiltonians -ff(O) and H{tf) are given, 
and set the initial and final configurations of the external parameters. Then we 
define I{t) and its eigenvectors accordingly, so that the commutation relations 
are obeyed at the boundary times and, finally, H{t) is designed via Eq. (6). 
While the an{t) may be taken as fully free time-dependent phases in principle, 
they may also be constrained by a pre-imposed or assumed structure of H{t). 
Sees. 3, 4 and 5 present examples of how this works for expansions, transport 
and internal state control. 

A generalization of this inverse method for non-Hermitian Hamiltonians was 
considered in Ibanez et al. (2011). Inverse engineering was applied to accel- 
erate the slow expansion of a classical particle in a time-dependent harmonic 
oscillator without final excitation. This system may be treated formally as a 
quantum two-level system with non-Hermitian Hamiltonian (Gao et al., 1991, 
1992). 

Quadratic in momentum invariants. - Lewis and Riesenfeld (1969) paid spe- 
cial attention to the time-dependent harmonic oscillator and its invariants 
quadratic in position and momentum. Later on Lewis and Leach (1982) found, 
in the framework of classical mechanics, the general form of the Hamiltonian 
compatible with quadratic-in-momentum invariants, which includes non har- 
monic potentials. This work, and the corresponding quantum results of Dhara 
and Lawande (1984) constitutes the basis of this subsection. 

A one-dimensional Hamiltonian with a quadratic-in-momentum invariant must 
have the form H = p^/2m + V{q,t),^ with the potential (Lewis and Leach, 

^ q and p may denote operators or numbers. The context should clarify their exact 
meaning. 
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1982; Dhara and Lawande, 1984) 



T77 1 

V{q,t) = -F{t)q + -u\t)q' + -^^U 



Pit) 



(7) 



p, Qc, u, and F are arbitrary functions of time that satisfy the auxihary equa- 
tions 



F{t)/m, 



(9) 



where ooq is a constant. Their physical interpretation will be explained below 
and depends on the operation. A quadratic-in-p dynamical invariant is given, 
up to a constant factor, by 

I = - mqc) - mp{q - qc)]'^ + ^muo (^^-^^ + U (^^-y^^ . (10) 

Now On in Eq. (4) satisfies (Lewis and Riesenfeld, 1969; Dhara and Lawande, 
1984) 



a„. = --l dt — + - 



h 



.p 



2p2 



:iii 



and the function </>„ can be written as (Dhara and Lawande, 1984) 

0„(g,t) = e^hV2p+(..P-..Pk/p]_L$„|^i^ (12) 



in terms of the solution $„(cr) (normalized in cr-space) of the auxiliary Schrodinger 
equation 



2m da"^ 2 



(13) 



The strategy of invariant-based inverse engineering here is to design p and q^ 
first so that I and H commute at initial and final times, except for launching 
or stopping atoms as in Torrontegui et al. (2011). Then H is deduced from 
Eq. (7). Applications will be discussed in Sees. 3 and 4. 
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2.2 Counterdiahatic or Transitionless Tracking Approach 



For the transitionless driving or counterdiahatic approach as formulated by 
Berry (2009), and equivalently by Demirplak and Rice (2003, 2005, 2008) ^ 
the starting point is a time-dependent reference Hamiltonian, 

Ho{t)=Y.\no{t))E^\t){n,{t)\. (14) 

n 

The approximate time-dependent adiabatic solution of the dynamics with Hq 
takes the form 

l^rW) = e^^"^*Vo(t)), (15) 
where the adiabatic phase reads 

Ut) = -\ f dt' E^^\t') + I fdt'inoit'MMt'))- (16) 
n Jo Jo 

The approximate adiabatic vectors in Eq. (15) are defined differently from 
the dynamical modes of the previous subsection, but they may potentially 
coincide, as we shall see. Defining now the unitary operator 

f/ = ^e*?"W|no(t))(rio(0)|, (17) 

n 

a Hamiltonian H{t) = ihUW can be constructed to drive the system exactly 
along the adiabatic paths of Ho{t), as H{t) = HQ{t) + Hcdit), where 



Hcdit) = ihY^ \\dMt)){no{t)\ - {nQ{t)\dMt))\no{t)){n,{t)\^ (18) 
is purely non-diagonal in the {\nQ{t))} basis. 

We may change the E^\t), and therefore -ffo(^) itself, keeping the same |no(t)). 
We could for example make all the E^\t) zero, or set ^„(t) = (Berry, 2009). 
Taking into account this freedom the Hamiltonian for transitionless driving 
can be generally written as 



H{t) = -hY,\nQ{t))Un^{t)\+ihY,\dtn^{t)){nQ{t)\. 

n n 

Subtracting Hcd{t), the generic Hq is 

Ho{t) = E \no{t))Uh{no{t)\dtno{t)) - hQ {no{t)\. 



(19) 



(20) 



^ Berry's transitionless driving method is equivalent to the counterdiahatic ap- 
proach of Demirplak and Rice (2003, 2005, 2008). In Section 2.4 we shall see how 
to further exploit this scheme together with "superadiabatic iterations" . 
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It is usually required that Hcdif) vanish for t < and t > tf, either suddenly 
or continuously at the boundary times. In that case the {\nQ{t))} become also 
at the extreme times (at least at t = 0~ and t = ) eigenstates of the full 
Hamiltonian. 

Using Eq. (1) and the orthonormality of the {|'^o(0))} we may write invariants 
of H(t) with the form I{t) = J2n \no(t))Xn{no(t)\. For the simple choice A„ = 
£;W(0), then J(0) = Ho{0). 

In this part and also in Sec. 2.1 the invariant-based and transitionless-tracking- 
algorithm approaches have been presented in a common language to make their 
relations obvious. Reinterpreting the phases of Berry's method as $,n{t) = 
and the states as |no(t)) = the Hamiltonians H{t) in Eqs. 

(6) and (19) may be equated. As well, the -ffo(^) implicit in the invariant- 
based method is given by Eq. (20), so that the dynamical modes can be also 
understood as approximate adiabatic modes of i^o(^) (Chen et al., 2011a). 
An important caveat is that the two methods could coincide but they do not 
have to. Given H{0) and H{tf), there is much freedom to interpolate them 
using different invariants, phase functions, and reference Hamiltonians HQ{t). 
In other words, these methods do not provide a unique shortcut but families 
of them. This flexibility enables us to optimize the path according to physical 
criteria and/or operational constraints. 

Non-Hermitian Hamiltonians. - A generalization is possible for non-Hermitian 
Hamiltonians in a weak non-hermiticity regime (Ibafiez et al., 2011, 2012b). It 
was applied to engineer a shortcut laser interaction and accelerate the decay 
of a two-level atom with spontaneous decay. Note that the concept of "popu- 
lation" is problematic for non-Hermitian Hamiltonians (Leclerc et al., 2012). 
This affects in particular the definition of "adiabaticity" and of the shortcut 
concept. It is useful to rely instead on normalization-independent quantities, 
such as the norm of a wave-function component in a biorthogonal basis (Ibanez 
et al., 2012b). 

Many-body Systems.- Following del Campo et al. (2012), the transitionless 
quantum driving can be extended as well to many-body quantum critical sys- 
tems, exploiting recent advances in the simulation of coherent A;-body inter- 
actions (Miiller et al., 2011; Barreiro et al., 2011). In this context STA allow 
a finite-rate crossing of a second order quantum phase transition without cre- 
ating excitations. Consider the family of quasi-free fermion Hamiltonians in 
dimension _D, 1-Lq = Z^kV^k Pk(A(t)) ■ a^] ipk, where the k-mode Pauli matrices 
are = (c"k' '^k' '^k) ^^"^ ^k — ('^ki''^k2) fermionic operators, and the 
sum goes over independent k-modes. Particular instances of quantum critical 
models within this family of Hamiltonians are the Ising and XY models in 
D = 1 (Sachdev, 1999), and the Kitaev model in D = 2 (Lee et al., 2007) and 
D = 1 (Sengupta et al., 2008). The function ak(A) = (O'k(A), a^(A), ak(A)) 
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is specific for eacli model (Dziarmaga et al., 2010). All these models can 
be written down as a sum of independent Landau-Zener crossings, where 
the instantaneous k-mode eigenstates have eigenenergies £k,± = il^fkl-^)! = 
±^Ja^{Xy + a^(A)2 + a^(A)2. It is possible to adiabatically cross the quan- 
tum critical point driving the dynamics along the instantaneous eigenmodes 
of Ho provided that the dynamics is driven by the modified Hamiltonian 
H = Ho + Tied, where (del Campo et al., 2012) 



Hcd = A'(t) E MM2 V'I [(«k(A) X dxM^)) ■ <Tk] (21) 

is typically highly non-local in real spaces and involves many-body interactions 
in the spin representation. However, it was shown in the ID quantum Ising 
model that a truncation of Tied with interactions restricted to range M is 
efficient to suppress excitations on modes k > (del Campo et al., 2012). 



2.3 Fast-forward Approach 



Based on some earlier results (Masuda and Nakamura, 2008), the fast-forward 
(FF) formalism for adiabatic dynamics and application examples were worked 
out in Masuda and Nakamura (2010, 2011); Masuda (2012) for the Gross- 
Pitaevskii equation or the corresponding Schrodinger equation. The aim of 
the method is to accelerate a "standard" system subjected to a slow varia- 
tion of external parameters by canceling a divergence due to an infinitely-large 
magnification factor with the infinitesimal slowness due to adiabaticity. A fast- 
forward potential is constructed which leads to the speeded- up evolution but, 
as a consequence of the different steps and functions introduced, the method is 
somewhat involved, which possibly hinders a broader application. The stream- 
lined construction of fast-forward potentials presented in Torrontegui et al. 
(2012a) is followed here. 

The starting point is the 3D time-dependent Gross-Pitaevskii (CP) equation 
(Dalfovo et al., 1999) 

= -|^^'^^^'^) + ^(x,t)^(x,t) +^3|^(x,t)|V(x,t). (22) 

Using the ansatz ■{/'(x, t) = r(x, t)e^'^^^'^'^ ('"(x, t), 0(x, t) G M) we formally solve 
for y(x, t) in (22) and get for the real and imaginary parts 
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Re[V{^, t)] = -h<P+^(^^- {V<pf ) - g,T\ (23) 



Im X, t = + — — + 

r 2m y r 

Imposing Im[V^(x, t)] = 0, i.e., 



(24) 



Eq. (23) gives a real potential. In the inversion protocol it is assumed that the 
full Hamiltonian and the corresponding eigenstates are known at the boundary 
times. Then we design r(x, t), solve for in Eq. (25), and finally get the 
potential V from Eq. (23). In Torrontegui et al. (2012a) it was shown how the 
work of Masuda and Nakamura (2008, 2010, 2011) relates to this streamlined 
construction. 



Since the phase </> that solves Eq. (25) depends in general on the particular 
r(x, t), Eq. (23) gives in principle a state-dependent potential. However, in 
some special circumstances, the fast-forward potential remains the same for 
all modes. This happens in particular for the Schrodinger equation, = 0, and 
Lewis-Leach potentials associated with quadratic-in-momentum invariants. In 
other words, the invariant-based approach can be formulated as a special case 
of the simple inverse method (Torrontegui et al., 2012a). 



2.4 Alternative Shortcuts Through Unitary Transformations 



Shortcuts found via the methods described so far or by any other approach 
might be difficult to implement in practice. In the cd approach, for instance, 
the structure of the complementary Hamiltonian Hcd could be quite different 
from the structure of the reference Hamiltonian Hq. Here are three examples, 
the first two for a particle of mass m in ID, the third one for a two-level 
system: 

- Example 1: Harmonic oscillator expansions (Muga et al., 2010), see Sec. 3: 
Ho = pV(2m) + mu\y2, = -{pq + qp)oo/i^uj). (26) 



- Example 2: Harmonic transport with a trap of constant frequency Uo/^n and 
displacement qo(t) (Torrontegui et al., 2011), see Sec. 4: 

Ho=py{2m) + {q-qo{t)fmu^j2, H^d = pqo- (27) 
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- Example 3: Population inversion in a two- level system (Berry, 2009; Chen 
et al., 2010c; Ibafiez et al., 2012a), see Sec. 5: 



Ho 








) 



(28) 



where 0o = arccos(Zo/-Ro) is a polar angle and Rq = {Xq + ZqY^'^. 

In all these examples the experimental implementation of Hq is possible, but 
the realization of the counter-diabatic terms is problematic. A way out is pro- 
vided by unitary transformations that generate alternative shortcut protocols 
without the undesired terms in the Hamiltonian (Ibanez et al., 2012a). A stan- 
dard tool is the use of different interaction pictures for describing one physical 
setting. Unitary operators W(t) connect the different pictures and the goal 
is frequently to work in a picture that facilitates the mathematical manipu- 
lations. In this standard scenario all pictures describe the same physics, the 
same physical experiments and manipulations. 

The main idea in Ibanez et al. (2012a) is to regard instead the unitary trans- 
formations as a way to generate different physical settings and different ex- 
periments, not just as mathematical transformations. The starting point is 
a shortcut described by the Schrodinger equation ihdtip(t) = H(t)ip(t), our 
reference protocol. (In all the above examples H = Hq + H^i-) The new 
dynamics is given by ihdtip'{t) = H' {t)ip' (t) , where ip'(t) = U{t)^ip{t), and 
H' = W{H - K)U, where K = ihUW. If U{0) = U{tf) = 1 the final states 
will coincide, i.e., ip'{tf) = ip(tf) for a given initial state ip'{0) = ip{0). If, in 
addition, W(0) = U{tf) = 0, then H{0) = H'{0), and H{tf) = H'{tf). Let 
us now list the unitary transformations that provide for the three examples 
realizable Hamiltonians (Ibanez et al., 2012a): 

- Example 1: Harmonic oscillator expansions, 



- Example 2: Harmonic transport, 

U = exp{-imqoq/h), H' = / {2m) + [q - %{t)fmujl/2, (30) 
where q'^ = qo + qo/ujl- 




(29) 



where a;' = _ g ^ ^ 



1/2 
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u = 








^ e*'^/^ ) 





Example 3: Population inversion in a two-level system, 

Zo - P 

.1, (31) 
P -Zq + h(l)/2 

where = arctan (^160/2X0), < < 27r, and P = [X^ + (;ieo/2)2] V^- 

Why do the Ws in Eqs. (29-31) have these forms? The answer lies in the sym- 
metry possesed by the Hamiltonian. Transformations of the form Li = e*-^^*^*^^ 
based on generators Gj of the corresponding Lie algebra produce operators 
within the algebra and, by suitably manipulating the function f{t) undesired 
terms may be eliminated. 

Superadiahatic iterations.- As discussed in Sec. 2.2, Demirplak, Rice, and 
Berry proposed to add a suitable counterdiabatic (cd) term^ Hf^ to the 
time dependent Hamiltonian -ffo(^) so as to follow the adiabatic dynamics 
of Hq. The same H^^ also appears naturally when studying the adiabatic 
approximation of the original system, i.e., the one evolving with Hq. This sys- 
tem behaves adiabatically, following the eigenstates of Hq., precisely when the 
counterdiabatic term is negligible. 

This is evident in an interaction picture (IP) based on the unitary transfor- 
mation Ao(t) = l^o(^))(?^o(0)| such that ipiit) = AqiPq. In this IP, the new 
Hamiltonian is Hi{t) = Al{t){Ho{t) - Ko{t))Ao{t) and Ko{t) = ihAo{t)Al{t). 
If Koit) is zero or negligible, Hi(t) becomes diagonal in the basis {|%(0))}, 
so that the IP equation is an uncoupled system with solutions 

mt)) = E |no(0))e-^/o ^"°'(*')^*'(no(0)|^i(0)). (32) 

n 

Correspondingly, \Mt)) = En |no(t))e-^/o ^""'(*')'^*'(no(0)|^o(0)). 

The same solution, which, for a non-zero Kq, is only approximate, is found ex- 
actly by adding to the IP Hamiltonian the counterdiabatic term ylQ(t) J<'o(t)Ao(t). 
This requires an external intervention and changes the physics of the original 
system. In the IP the modified Hamiltonian is H^^^ = Hi + Al{t)Ko{t)Ao{t) = 
Ali{t)Ho{t)Ao{t) and in the Schrodinger picture (SP) the modified Hamilto- 
nian IS H^^\t) = Ho{t) + Ko{t), so we identify H^^it) = Ko{t). In other words, 
a "small" coupling term Kq that makes the adiabatic approximation a good 
one also implies a small counterdiabatic manipulation. However, irrespective 
of the size of Kq, HQ^\t) provides a shortcut to slow adiabatic following be- 
cause it keeps the populations in the instantaneous basis of Hq invariant, in 



^ This is the H^d term of Sec. 2.2. The superscript is added now to distinguish 
it from higher order cd terms introduced below. 
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particular at the final time tf. 

Looking for generalized adiabatic approximations, Garrido (1964), Berry (1987) 
or Descliamps et al. (2008) have investigated further iterative interaction pic- 
tures and the corresponding approximations. The idea is best understood by 
working out explicitly the next iteration: one starts with ihdtipi{t) = Hi%l)i{t) 
and diagonalizes Hi{t) to produce its eigenbasis {|ni(t))}. A unitary operator 
^1 = |'^i(^))('^i(0)| plays now the same role as Aq in the previous IP. It 
defines a new IP wave function ip2{t) = ^i(^)V'i that satisfies ihdtip2it) = 
H2Mt), where H2it) = A\{t){Hi{t) - Ki{t))Ai{t) and A'l = zMi4- If 
is zero or "small" enough, i.e. if a (first order) superadiabatic approximation 
is valid, the dynamics would be uncoupled in the new interaction picture, 
namely, 

mt)) = E l^i(0))e-^/o^^''(*')'^*' (^^(0)1^.2(0)). (33) 

We may get the same result by changing the physics and adding A[(t)Ki(t)Ai(t) 
to H2 (Demirplak and Rice, 2008; Ibafiez et al., 2012a). In the SP the added 
interaction becomes a first order counterdiabatic term Hj.^^ = AqKiAq. Trans- 
forming back to the SP and using Aj{0) = 1 the state (33) becomes 

l^o(t)) = EE l"^o(t))(mo(0)|ni(t))e-^/o ^^''{*')'^*'(ni(0)|^o(0)). (34) 

71 m 

Quite generally the populations of the final state in the adiabatic basis {\no(tf))} 
will be different from the ones of the adiabatic process, unless |no(0)) = 
\ni(tf)) and |ni(0)) = |no(0)), up to phase factors. The first condition is 
satisfied if -K'o(^/) = and the second one if Kq{0) = 0. Then the su- 
peradiabatic process will actually lead to the same final populations as an 
adiabatic one, possibly with different phases for the individual components. 
Similarly, the first-order counterdiabatic term H^^ would provide a shortcut 
with Hq^^ = Hq + H^^ in the SP, different from the one carried out by Hq^\ 
Moreover, if Ki{0) = Ki{tf) = 0, then H^^^ = Hq, at t = 0,t/. Further 
iterations define higher order superadiabatic frames. Is there any advantage 
in using one or another counterdiabatic scheme? There are two reasons that 
could make higher order schemes attractive in practice: one is that the struc- 
ture of the H^'i^ may change with j. For example, for a two-level population 
inversion H^'^ = h{Qo/2)ay, whereas if^^ = ^(9i/2)(cos GqCx — sinOoCTz), 
where Gi is the polar angle corresponding to the the Cartesian components 
of Hi = Xiax + YiO'y + Zia^ (Ibanez et al., 2012a). The second reason is that, 
for a fixed process time, the cd-terms are smaller in norm as j increases, up 
to a value in which they begin to grow, see e.g. Deschamps et al. (2008). One 
should pay attention though not only to the size of the cd-terms but also to 
the feasibility of the boundary conditions at the time edges to really generate 
shortcuts in this manner. 
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2.5 Optimal Control Theory 



Optimal control theory (OCT) is a vast field covering many techniques and 
apphcations. As for STA, fast expansions (Salamon et al., 2009), wavepacket 
splitting (Hohenester et al., 2007; Grond et al., 2009a,b), transport (Murphy 
et al., 2009) and many-body state preparation (Rahmani and Chamon, 2011) 
have been addressed with different OCT approaches. The combination of OCT 
techniques with invariant-based engineering STA is particularly fruitful since 
the later provides by construction families of protocols that achieve a per- 
fect fidelity or vanishing final excitation, whereas OCT may help to select 
among the many possible protocols the ones that optimize some physically 
relevant variable (Stefanatos et al., 2010, 2011; Chen et al., 2011b; Stefanatos 
and Li, 2012). In this context the theory used so far is the maximum princi- 
ple of Pontryagin (1962). For a dynamical system x = f(x(t),n), where x is 
the state vector and u the scalar control, in order to minimize the cost func- 
tion J{u) = Jq^ g{'x{t),u)dt, the principle states that the coordinates of the 
extremal vector x(t) and of the corresponding adjoint state p{t) formed by 
Lagrange multipliers, fulfill Hamilton's equations for a control Hamiltonian 
He = Po5'(x(t), m) + p"^ ■ f (x(t), m). For almost all times during the process He 
attains its maximum at m = u{t) and He = c, where c is constant. We shall 
discuss specific applications in Sees. 3 and 4. 



3 Expansions of trapped particles 

Performing fast expansions of trapped cold atoms without losing or exciting 
them is important for many applications: for example to reduce velocity dis- 
persion and coUisional shifts in spectroscopy and atomic clocks, decrease the 
temperature, adjust the density to avoid three body losses, facilitate temper- 
ature and density measurements, or to change the size of the cloud for further 
manipulations. Of course trap compressions are also quite common. 

For harmonic traps we may address expansion or compression processes with 
the quadratic- in-p invariants theory by setting qc = U = F = in Eq. (7). This 
means that Eq. (9) does not play any role and the important auxiliary equation 
is the "Ermakov equation" (8). The physical meaning of p is determined by 
its proportionality to the standard deviation of the position of the "expanding 
(or contracting) modes" 

Here we shall discuss the expansion from u}{0) = Uq to uj{tf) = Uf (Chen et 
al., 2010b). Choosing 

p(0) = l, p(0) = 0, (35) 
H[0) and /(O) commute. They actually become equal, and have common eigen- 
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functions. Consistent with the Ermakov equation, p(0) = holds as well for a 
continuous frequency. At we impose^ 

p{t^)=^ = {uj,/ujjf'\ p(t/) = 0, p(t/) = 0. (36) 

In this manner the expanding mode is an instantaneous eigenvector of H at 
t = and tf, regardless of the exact form of pit). To fix p(t), one chooses 
a functional form to interpolate between these two times, flexible enough 
to satisfy the boundary conditions. For a simple polynomial ansatz p{t) = 
6 (7 - 1) (7 - 1) s-^+lO (7 - 1) s^+l (Palao et al., 1998), where s = t/tj. 

The next step is to solve for uj(t) in Eq. (8). This procedure poses no funda- 
mental lower limit to t/, which could be in principle arbitrarily small. There 
are nevertheless practical limitations and/or prices to pay. For short enough 
tf, u{t) may become purely imaginary at some t (Chen et al., 2010b) and the 
potential becomes a parabolic repeller. Another difficulty is that the transient 
energy required may be too high, as discussed in Chen and Muga (2010) and in 
the following subsection. Since actual traps are only approximately harmonic, 
large transient energies will imply perturbing effects of anharmonicities and 
thus undesired excitations of the final state, or even atom losses. 

3.1 Transient Energy Excitation 

Knowing the transient excitation energy is also important to quantify the prin- 
ciple of unattainability of zero temperature, first enunciated by Nernst. This 
principle is usually formulated as the impossibility to reduce the temperature 
of any system to the absolute zero in a finite number of operations, and identi- 
fied with the third law of thermodynamics. Kosloff and coworkers in (Salamon 
et al., 2009) have restated the unattainability principle in quantum refriger- 
ators as the vanishing of the cooling rate when the temperature of the cold 
bath approaches zero, and quantify it by the scaling law that relates cooling 
rate and cold bath temperature. We shall examine here the consequences of 
the transient energy excitation on the unattainability principle in two ways: 
for a single, isolated expansion, and considering the expansion as one of the 
branches of a quantum refrigerator cycle (Chen and Muga, 2010). 

A lower bound S„ for the time-averaged energy of the n-th expanding mode 
En (time averages from to will be denoted by a bar) is found by applying 

^ If p{tf) / the final frequency would not be ujf but uj{tf) = [uij- — p/'j]^^'^. If 
discontinuities are allowed and the frequency is changed abruptly from 00 [tj) to ujf 
the excitations will also be avoided, at least in principle. A similar discontinuity 
is possible at t = if /3(0) 7^ and the frequency jumps abruptly from ujq to 
a;(0) = [a;2-p(0)]V2. 
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calculus of variations (Chen and Muga, 2010), so that > S„. If the final 
frequency is small enough to satisfy t/ ^ I/^wqW/, and 7^1, the lower 
bound has the asymptotic form ~ {2n + l)h/{2ujft'j). A consequence is 

that tf > {2n + l)h/ {2ujfEn)- When En is limited, because of anharmonic- 
ities or a finite trap depth, the scaling is fundamentally the same as the one 
found for bang-bang methods with real frequencies (Salamon et al., 2009), 
and leads to a cooling rate R oc Tj^^'^ in an inverse quantum Otto cycle (the 
proportionality factor may be improved by increasing the allowed En). This 
dependence had been previously conjectured to be a universal one character- 
izing the unattainability principle for any cooling cycle (Rezek et al., 2009). 
The results in Chen and Muga (2010) provide strong support for the validity 
of this conjecture within the set of processes defined by ordinary harmonic os- 
cillators with time-dependent frequencies. In (Hoffmann et al., 2011) a faster 
rate ~ —Tc/ log Tc is found with optimal control techniques for bounded trap 
frequencies, allowed to become imaginary. There is no contradiction with the 
previous scaling since bounding the trap frequencies does not bound the sys- 
tem energy. In other words, achieving such fast cooling is not possible if the 
energy cannot become arbitrarily large. 

Independently of the participation of the harmonic trap expansion as a branch 
in a refrigerator cycle, we may apply the previous analysis also to a sin- 
gle expansion, assuming that the initial and final states are canonical den- 
sity operators characterized by temperatures Tq and Tj. These are related by 
Tf = {u!f/uJo)To for a population-preserving process. In a harmonic potential 
expansion, the unattainability of a zero temperature can be thus reformulated 
as follows: The transient excitation energy becomes infinite for any population- 
preserving and finite-time process when the final temperature is zero (which 
requires Uf = 0). The excitation energy has to be provided by an external 
device, so a fundamental obstruction to reach Tj = in a finite time, is the 
need for a source of infinite power (Chen and Muga, 2010). 

The standard deviation of the energy was also studied numerically (Chen and 
Muga, 2010). There it was found that the dominant dependences of the time 
averages scale with Uf and t/ in the same way as the average energy. These 
dependences are different from the ones in the Anadan and Aharonov (1990) 
(AA) relation AHtf > |, where AH = /q^ AH{t)dt/tf. 

3.2 Three Dimensional Effects 

The previous discussion is limited to one dimension (ID) but actual traps 
are three-dimensional and at most effectively ID. Torrontegui et al. (2012c) 
worked out the theory and performed numerical simulations of fast expan- 
sions of cold atoms in a three-dimensional Gaussian-beam optical trap. Three 
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different methods to avoid final motional excitation were compared: inverse 
engineering using Lewis-Riesenfeld invariants, which provides the best overall 
performance, a bang-bang approach with one intermediate frequency, and a 
"fast adiabatic approach" ^ . 

The optical trap considered in Torrontegui et al. (2012c) is formed by a laser, 
red detuned with respect to an atomic transition, and is characterized in the 
harmonic approximation by longitudinal and radial frequencies. To fourth or- 
der in the coordinates the effective potential includes anharmonic terms and 
radial-longitudinal coupling terms. While magnetic traps allow for an inde- 
pendent control of longitudinal and radial frequencies (Schaff et al., 2010, 
2011a,b), this is not the case for a simple laser trap. In Torrontegui et al. 
(2012c) it was assumed that the time- dependence of the longitudinal frequency 
is engineered to avoid final excitations with a simple ID harmonic theory. The 
main conclusion of the study is that the transitionless expansions in opti- 
cal traps are feasible under realistic conditions. For the inverse engineering 
method, the main perturbation is due to the possible adiabaticity failure in 
the radial direction, which can be suppressed or mitigated by increasing the 
laser waist. This waist increase would also reduce smaller perturbing effects 
due to longitudinal anharmonicity or radial-longitudinal coupling. The simple 
bang-bang approach fails because the time for the radial expansion is badly 
mismatched with respect to the ideal time, and the fast adiabatic method 
fails for short expansion times as a result of longitudinal excitations. Compli- 
cations such as perturbations due to different noise types, and consideration of 
condensates, gravity effects, or the transient realization of imaginary trap fre- 
quencies are still open questions. Other extensions of Torrontegui et al. (2012c) 
could involve the addition of a second laser for further control of the potential 
shape, or alternative trap shapes. Optical traps based on Bessel laser beams, 
for example, may be useful to decouple longitudinal and radial motions. 



3.3 Bose-Einstein Condensates 

In this section we shall discuss the possibility of realizing STA in a harmon- 
ically trapped Bose-Einstein condensate using a scaling ansatz. A mean-field 
description of this state of matter is based on the time-dependent Gross- 
Pitaevskii equation (GPE) (Dalfovo et al., 1999), 



The adiabaticity condition for the harmonic oscillator is |\/2w/(8a;^)| <^ 1. An 
efficient, but still adiabatic, strategy by Chen et al. (2010b) is to distribute w/w^ 
uniformly along the trajectory, i.e., oj/oo'^ = c, c being constant. Solving this differ- 
ential equation and imposing uj = oj{tf) we get a;(t) = ujq/[1 — (uj — uJo)t/{tf(jjf)]. 
This may be enough for some applications. This function was successfully applied 
in Bowler et al. (2012). 
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9^(x,t) 
dt 



1 

—A + -mu^it)^^ + (7z3|^(x,t)p 



^(x,t). (37) 



Here, A is the D- dimensional Laplacian operator and (^d is the D- dimensional 
coupling constant. For a three-dimensional cloud, using the normalization 
/ |^(x, i)p(ix = 1, (73 = ^^^^^ , for a condensate of a number of atoms 
of mass m, interacting with each other through a contact Fermi-Huang pseu- 
dopotential parameterized by a s-wave scattering length a. In D = 1,2 the 
corresponding expression for can be obtained by a dimensional reduction of 
the 3D GPE (Salasnich et al., 2002). As a mean- field theory the GPE overesti- 
mates the phase coherence of real Bose-Einstein condensates. The presence of 
phase fluctuations generally induces a breakdown of the dynamical self-similar 
scaling law that governs the dynamics of the expanding cloud and the forma- 
tion of density ripples. The conditions for quantum phase fluctuations to be 
negligible for STA were discussed in del Campo (2011a). In the following we 
shall ignore phase fluctuations. The results of this section will be generalized 
to strongly correlated gases in Sec. 3.4, including as a particular case, the 
microscopic model of ultracold bosons interacting through s-wave scattering. 

STA in the mean-field regime were designed in Muga et al. (2009) based 
on the classic results by Castin and Dum (1996), and Kagan, Surkov and 
Shlyapnikov (1996), who found the exact dynamics of the condensate wave- 
function under a time-modulation of the harmonic trap frequency. Consider 
a condensate wavefunction \E'(x, t = 0), a solution of the time-independent 
GPE with chemical potential /i in a harmonic trap of frequency uq, i.e., 
(-^A + imw^x^ + gol^i^^t = 0)\^ - /i)^(x,t = 0) = 0. Under a mod- 
ulation of the trap frequency u{t) the scaling ansatz 



^(x,t) 



1 

— exp 

p2 



.m|x| p 
^ 2h ~p 



./ir(t) 

' h 



,t = 



(3^ 



is an exact solution of the time-dependent Gross-Pitaevskii equation provided 
that 



p + c.(i)V = 4. = r(t)=f%. (39) 

p-^ p"^ ^ Jo p 

It follows that the scaling factor p must be a solution of the Ermakov equa- 
tion, precisely as in the single-particle harmonic oscillator case. This paves 
the way to engineer a shortcut to an adiabatic expansion or compression from 
the initial state \E'(x, t = 0) to a target state \E^(x, tj) = \E^(x/p, t = 0)/p~ 
by designing the trajectory p(t). The modulation of the coupling constant 
required in D = 1,3 can be implemented with the aid of a Feshbach reso- 
nance (Muga et al., 2009), or, in D = 1, by a modulation of the transverse 
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confinement (Staliunas et al., 2004; Engels et al., 2007; del Campo, 2011a). 
The D = 2 requires no tuning in time of tlie coupling constant as a result of 
the Pitaevskii-Rosch symmetry (Pitaevskii and Rosch, 1997). It has recently 
been suggested that this symmetry is broken upon quantization, constituting 
an instance of a quantum anomaly in ultracold gases (Olshanii et al., 2010). 
To date no experiment has provided evidence in favor of this observation. We 
point out that observing a breakdown of shortcuts to expansions of 2D BEC 
clouds would help to verify this quantum-mechanical symmetry breaking. 

An important simplification occurs in the Thomas-Fermi regime, where the 
mean-field energy dominates over the kinetic part. Assuming the validity of 
this regime along the dynamics, the scaling ansatz (38) becomes exact as long 
as the following consistency equations are satisfied. 



* dt' 



p + uj{tYp = -j^, gn{t) =gn{t = 0), r(t) = / — . (40) 
p^+^ Jo p^ 

Hence, in the Thomas-Fermi regime, it is possible to engineer a shortcut ex- 
actly, while keeping the coupling strength gn constant (Muga et al., 2009). 
Optimal control theory has been recently applied in this regime to find opti- 
mal protocols with a restriction on the allowed frequencies Stefanatos and Li 
(2012). 

Dimensional reduction and modulation of the non-linear interactions.- For 
low dimensional BECs, tightly confined in one or two directions, an effective 
tuning of the coupling constant can be achieved by modulating the trapping 
potential along the tightly confined axis, see e.g. Staliunas et al. (2004), a 
proposal experimentally explored in Engels et al. (2007). In a nutshell, the 
tightly confined degrees of freedom decoupled from the weakly confined ones, 
are governed to a good approximation by a non-interacting Hamiltonian. It 
is then possible to perform a dimensional reduction of the 3D GPE, and de- 
rive a lower- dimensional version for the weakly confined degrees of freedom, 
where the effective coupling constant inherits a dependence of the width of 
the transverse modes which have been integrated out. Adiabatically tuning 
the transverse confinement leads to a controlled tuning of the effective cou- 
pling constant. A faster-than-adiabatic modulation can be engineered by im- 
plementing a shortcut in the transverse degree of freedom. Consider the 3D 
mean-field description 



ih 



a^(x,t) 



2 



;^A + V<=^(x,t) + (73|*(x,t)P *(x,t), (41) 
2m J 



dt 

with V'^''(x, t) = '^[ux{tYx^ + Uy{t)'^y'^+Uz{tYz'^]. For tight transverse confine- 



ment {ojx ^ iOy ^ ujz and N\a\JmuJz/h ^ 1), the transverse excitations are 
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frozen. The transverse mode can be approximated by the single-particle har- 
monic oscillator ground state ^Q{x,y,t), so that the wavefunction factorizes 
\l/(x, t) = ^Q{x,y,t)ip{z,t). Integrating out the transverse modes, and up to a 
time-dependent constant which can be gauged away, one obtains the reduced 
GPE 



ih 



dt 



2m dz"^ 



(42) 



with the effective coupling gi(t) = JJdxdy\^o{x,y,t)\^. A general trajectory 
gi{t) can be implemented by modifying the frequency u}±{t) of the transverse 
confinement according to 



u;l{t)=ul{0) 



I hit) 
2giit) 



9iit) 



(43) 



in quasi-lD atomic clouds (del Campo, 2011a). The first term in the RHS cor- 
responds to the adiabatic tuning discussed in Staliunas et al. (2004); Engels 
et al. (2007) while the remaining terms are associated with the STA dynamics 
in the transverse modes. A similar analysis applies to the control of the effec- 
tive coupling constant in a pancake condensate, in the x-y plane, under tight 
confinement along the ^-direction (del Campo, 2011a). 

We note that this technique is restricted to tune the amplitude of the cou- 
pling constant, at variance with alternative techniques based on Feschbach 
or confinement-induced resonances which can change both the amplitude and 
character of the interactions, e.g., from attractive to repulsive (Bloch et al., 
2008). 



3.4 Strongly Correlated Gases 



The preceding sections were focused on single-particle systems and a mean- 
field description of Bose- Einstein condensates. We have seen that the inversion 
of scaling laws is a powerful technique to design STA in those processes where 
the dynamics is self-similar, e.g. expansions, or transport. In the following 
we focus on the engineering of STA in strongly correlated quantum fluids of 
relevance to ultracold gases experiments. We shall consider a fairly general 
model in dimension D consisting of indistinguishable particles with coor- 
dinates Xj G M^, trapped in a time- dependent isotropic harmonic potential of 
frequency uj(t) and interacting with each other through a two-body potential 
V(xj — Xj). The many-body Hamiltonian describing this system reads (del 
Campo, 2011b) 
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N 
i=l 



2m 



■A,; 



-ma;^(t)x- 



+ e^V(xi -Xj), 



(44) 



where Aj is the D- dimensional Laplacian operator for the Xj variable, and 
e = e{t) is a dimensionless time-dependent coupling strength satisfying e(0) = 
1. We shall further assume that V(Ax) = A~"V(x) under scaling of the co- 
ordinates. Specific realizations of this model include the Calogero-Sutherland 
model (Sutherland, 1998), the Tonks- Girardeau gas (Ohberg and Santos, 2002; 
Minguzzi and Gangardt, 2005), Lieb-Liniger gas (Bujan et al., 2008), Bose- 
Einstein condensates (BEG) (Gastin and Dum, 1996; Kagan, Surkov and 
Shlyapnikov, 1996; Muga et al., 2009), including dipolar interactions (O'Dell 
et al., 2004), and more general many-body quantum systems (Gritsev et al., 
2010). For simplicity, we leave out other cases to which similar techniques can 
be applied, such as strongly interacting mixtures (Minguzzi and Girardeau, 
2007) or systems with internal structure (Mousavi et al., 2007; Deuretzbacher 
et al., 2008). 

Let us now consider an equilibrium state $ of the system (44) at t = 
with chemical potential /i. For compactness we shall use the notation ^j-k = 
{xj,Xj+i, . . . , Xk-i, Xk}. It is possible to find a self-similar scaling solution of 
the form 



$(Xl:iV,t) 



P 



D/2 



exp 



N 1 • 



E 



{ 2ph 



ifiT{t)/h 



$(^,t = 0j, (45) 



where r(t) = Jq dt' / p'^{t'), whenever the scaling factor p = p(t) is the solution 
of the Ermakov differential equation, p + u'^(t)p = ouq/p^, with ujq = cj(0), 
satisfying the boundary conditions p(0) = 1 and p(0) = 0. This is the same 
consistency equation that arises in the context of the single-particle time- 
dependent harmonic oscillator. 

Scaling laws greatly simplify the dynamics of quantum correlations. Let us 
consider the time-evolution of the n-particle reduced density matrix 



^„(xi,^;x;^„;t) = _^y j [l dXi^*{xi.,N;t)^{x.[.^,y^n+i:N;t). (46) 



i=n+l 



Provided the scaling law holds, its time evolution reads 



^„(x,..;x;^„;t) = p-V[ ^; ^; ) exp(-lA ^"=i(^; ) \ ^^^^^ 



P P 



P Wo 



where /q = ■Jh/mujQ. 
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Local quantum correlations depend exclusively on the diagonal elements of 
5'n(xi:n; ^i-n'i t) manifest directly the self-similar dynamics. For instance, 
the time evolution of the density profile n(x) = (7i(x;x) reads n(x, t) = 
p-"^n(f,t = 0). 

The dynamics of non-local correlations is more involved due to the presence 
of the oscillatory phase. As an example, the evolution of the one-body re- 
duced density matrix (OBRDM) under self-similar dynamics (Minguzzi and 
Gangardt, 2005; Gritsev et al, 2010), 

1 /xy \ /ipx^ — 

^i(x,y;t) = — ^1 -,-;0 exp ^72~ ' (^8) 

\P P J V P^o 2/^ / 

induces a non-self-similar evolution of the momentum distribution, its Fourier 
transform ?7.(k, t) = / dxdy e''^'^^~^^gi{x,y]t). It is expected that the oscil- 
latory phases distort quantum correlations. The case of a free expansion, 
where the frequency modulation in terms of the Heaviside function G(t) reads 
u{t) = UQQ{—t), has received much attention. The solution to the Ermakov 

equation for the scaling factor is p{t) = \Jl + io^t"^ and for t ^ cuq^, p{t) ~ ujot, 
p = Uq. Using the method of the stationary phase, it follows that 



n(k, t) ~ \27ruJolt/pr9i{uJo^lo/p, oookl'Jp), (49) 

i.e., the asymptotic momentum distribution is mapped to the scaled density 
profile of the initial state (Jukic et al., 2008; Oezer et al., 2009; Gritsev et al., 
2010). As a result, all information of the off-diagonal elements of the OBRDM 
is lost. Similar effects result in an expansion in finite-time tj ^ luq^ and signal 
the breakdown of adiabaticity. Excitations manifest as well in local correlation 
functions, e.g. excitation of the breathing mode of the cloud. 

In the adiabatic limit (r ^ (^o^), the time- variation of the the scaling factor 
vanishes p(t) ^ 0, resulting in the adiabatic trajectory p{t) = ujQ/uj{t). At 
all times the time-evolution of the OBRDM and the momentum distribution 
can be related by a scaling transformation of their form at t = 0, 

(?i(x, y; t) = -^g, (^-^, -|y; 0^ , n(k, t) = p^(t)n(p(t)k, 0). (50) 

These expressions can be applied for expansions (p(t) > 1) and compressions 
(p(t) < 1), and generally still require tuning the interaction coupling strength. 
Nonetheless, the required adiabatic time scale can be exceedingly long and we 
next tackle the problem of achieving a final scaled state in a predetermined 
expansion time tf. The upshot of the frictionless dynamics is that quantum 
correlations at the end of the quench {t = tf, and only then) are those of the 
initial state scaled by a factor p(t/) = 7 (del Campo, 2011b). In particular. 
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^7i(x, y-tf) = ( -, -; I , n{k,tf) = ^''nijk, 0). (51) 

7 \7 7 / 

Similar expressions hold for higher-order correlations, i.e. gn{'^i:n,yi:n]tf) = 
l^^Qn (xi:„/7, yi:n/7; 0). Moreover, as long as the initial state is an equilibrium 
state in the initial trap, so it is the state at tj with respect to the final trap, 
preventing any non-trivial dynamics after the quench, for t > tf if uj{t > tf) = 
ujf. Nonetheless, at intermediate times t G [0,t/) the momentum distribution 
exhibits a rich non-equilibrium dynamics, and can show for instance, evolution 
towards the scaled density profile of the initial state. 

We close this section with two comments. First, the applicability of STA based 
on inversion of scaling laws is not restricted to fermionic or bosonic systems, 
but can be applied as well to anyonic quantum fiuids for which dynamical 
scaling laws are known (del Campo, 2008). Systems with quantum statistics 
smoothly extrapolating between bosons and fermions might be realized in the 
laboratory following Keilmann et al. (2011). Second, the possibility of scaling 
up the system while preserving quantum correlations constitutes a new type 
of microscopy of quantum correlations in quantum fiuids (del Campo, 2011b; 
del Campo and Boshier, 2012). It is as well of interest to design new protocols 
to reconstruct the initial quantum state of the system from the time-evolution 
of its density profile (Bertrand and Betrand, 1987; Leonhardt and Schneider, 
1997), a tomographic technique demonstrated experimentally in Kurtsiefer et 
al. (1997), and apphcable to many-body systems (del Campo et al., 2008b). 

Scaling laws in other trapping potentials.- Scaling laws for many-body sys- 
tems can be found for more general types of confinements. Among them, ho- 
mogeneous potentials are of particular interest, since they simplify the corre- 
spondence between ultracold atom experiments and condensed matter theory. 
The early experimental implementations of the paradigmatic particle in a box 
aimed at the creation of optical billiards for ultracold gases (Milner et al., 2001; 
Friedman et al., 2001). Trapping of a BEC in an all-optical box was reported 
in Meyrath et al. (2005) and analogous traps have been created in atom chips 
(van Es et al., 2010). For the purpose of implementing STA, the dynamical 
optical dipole potential may be realized using the highly versatile "painting 
technique" , which creates a smooth and robust time-averaged potential with a 
rapidly- moving laser beam (Henderson et al., 2009), or alternatively, by spatial 
light modulators (Boyer et al., 2006). 

The breakdown of adiabaticity in a time-dependent homogeneous potential 
leads to quantum transients related to the diffraction in time (DIT) effect, 
see del Campo et al. (2009) for a review. A sharply localised matter-wave in 
a region of space, after sudden removal of the confinement, exhibits during 
free evolution density ripples. The earliest example discussed by Moshinsky 
(1952) the free evolution of a truncated cut-off plane wave, exhibits an oscil- 
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latory pattern with the same functional form than the diffraction pattern of 
a classical light beam from a semi-infinite plane. The phenomenon is ubiq- 
uitous in matter-wave dynamics induced by a quench, and in particular, it 
arises in time- dependent box potentials in one (Gerasimov and Kazarnovskii, 
1976; Godoy, 2002; del Campo and Muga, 2006), two and three (Godoy, 2003) 
dimensions. The effect manifests as well in strongly interacting gases such as 
ultracold bosons in the To nks- Girardeau regime (del Campo and Muga, 2006; 
del Campo, 2008). Moreover, when the piston walls move at a finite speed f , 
the adiabatic limit is not approached monotonically as f — )■ 0. It was shown in 
del Campo et al. (2008a); Mousavi (2012a,b) that an enhancement of DIT oc- 
curs when the walls move with the dominant velocity component of the initial 
confined state, due to a constructive interference between the expanding and 
reflected components from the walls. When the reflections from the confine- 
ment walls dominate, the non-adiabatic dynamics in time-dependent homoge- 
neous potentials lead to Talbot oscillations and weave a quantum carpet in the 
time evolution of the density profile (Friesch et al., 2000; Ruostekoski et al., 
2001). Suppression of these excitations is dictated by the adiabatic theorem 
both in the non-interacting (Berry and Klein, 1984; Dodonov et al., 1993; Chen 
et al., 2009; Mostafazadeh, 2001) and mean-field regime (Band et al., 2002). 
Further, for non-interacting systems one can prove that no shortcut based on 
invariants or scaling exists in time- dependent homogeneous potentials. At the 
single-particle level, this follows from the fact that the family of trajectories 
for the width ^ (t) of a box-like potential for which a dynamical invariant exist 
(Berry and Klein, 1984), takes the form ^(t) = [at^ 4- fet -|- cja, which is in- 
compatible with the boundary conditions required to reduce a time-evolving 
scaling solution to the initial and target states. For many-body quantum flu- 
ids, the same result is derived from the consistency equations for self-similar 
dynamics to occur. To find a shortcut in this scenario one has to relax the con- 
dition on the confinement and allow for a inhomogeneous auxiliary harmonic 
potential of the form (del Campo and Boshier, 2012) 



where x G M'^, |x| G [0,.^(t)]. For D = 1, a box with one stationary wall 
at a; = and moving wall at x = ^{t) is assumed. Cylindrical and spherical 
symmetry is imposed for D = 2, 3 respectively. This auxiliary potential can 
be implemented by means of a blue-detuned laser (Khaykovich et al., 2002) or 
direct painting with a rapidly moving laser (Milner et al., 2001; Friedman et 
al., 2001; Henderson et al., 2009). Thanks to its presence it is possible to find 
dynamical self-similar solutions to single-particle and many-body Schrodinger 
equations for time-dependent box-like confinements with a general modulation 
of the width ^{t). In particular, consider the Hamiltonian 




(52) 
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^ r ^2 

^ =E - + +eE V^(x. - X,), (53) 



i=l 



where Xj G M'^, = |x|j G [0,.^(t)], and let us introduce the scahng fac- 
tor p{t) = ^(t)/{(0). If V(Ax) = A-°V(x), e{t) = p(t)"-^ in the presence 
of t/^^^(x, t), the time evolution of an initial eigenstate of the system with 
chemical potential /i follows a scaling law in Eq. (45). Given the existence 
of a scaling law, a many-body shortcut can be engineered by designing the 
scaling factor as for the simple harmonic oscillator, ensuring that the time- 
evolving state reduces to the initial and target state at the begining and end 
of the evolution (del Campo and Boshier, 2012). Naturally, this is possible 
as well for Bose-Einstein condensates in the mean-field regime, extending 
the Castin-Dum-Kagan-Surkov-Shlyapnikov scaling ansatz (del Campo and 
Boshier, 2012). 

Along a shortcut to an adiabatic expansion, the auxiliary potential is expul- 
sive in an early stage of the expansion, expelling the atoms from the center 
and providing the required speed-up. The rapidly expanding cloud is slowed 
down in a second state of the expansion, when t/^"^(x, t) becomes a trapping 
potential. The sequence is reversed in a shortcut to an adiabatic compression. 
In both cases, at t = tf, f/'^"^(x, t) vanishes, and the cloud reaches the target 
state, a stationary state of the final Hamiltonian. As a result, STA provide a 
variant of the paradigmatic model of a quantum piston (Quan and Jarzynski, 
2012). 



3.5 Experimental Realization 



Experiments of fast shortcut expansions have been realized at Nice with mag- 
netic confining of ^''Rb atoms for ultra-cold clouds (Schaff et al., 2010) and 
condensates in the Tomas-Fermi regime (Schaff et al., 2011a). Compared to 
the simple expansions treated in (Chen et al., 2010b), gravity introduces and 
extra linear term in the Hamiltonian and requires a treatment with additional 
boundary conditions. 

For the cold cloud, samples of = 10^ atoms and temperature Tq = 1.63 
/iK were used to keep the time between collisions small ~ 28 ms, and the 
potential effectively harmonic. The initial trap frequencies for x, y, z directions 
in Hz were (228.1,22.2,235.8) and the final ones (18.1,7.1,15.7). The results 
for the fast (35 ms) 15-fold frequency decompression to the trap in the vertical 
dimension, yielded a residual center-of-mass oscillation of the cloud equivalent 
to that of a 1.3-s-long linear decompression, a reduction by a factor of 37. 

For the condensate, the number of atoms was A^ = 1.3 x 10^ and the ini- 
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tial temperature was Tq = 130 nK (Schaff et al., 2011a). The potential is 
U{r,t) = mu\(t){x'^ + z^) + |ma;||(t)t/^ + mgz. Initial radial {x,z) and ax- 
ial (y) frequencies were 235.8 Hz and 22.2 Hz, respectively. The experiment 
performed a 30-ms-long radial decompression of the trap by a factor of 9, 
yielding a final radial frequency of 26.2 Hz. The axial frequency was reduced 
by a factor of 3 to a final value 7.4 Hz. Using scaling techniques similar to the 
ones in Sec. 3.3 it was shown that this decompression is a shortcut for both 
directions. Residual excitations were attributed to imperfect implementation 
of uj{t), anharmonicities, and trap tilting. 

3.6 Optimal Control 

The time- dependent frequency of a harmonic trap expansion based on invari- 
ants can be optimized with respect to time or to transient excitation energy, 
restricting the allowed transient frequencies (Stefanatos et al., 2010, 2011). 

Kosloff and coworkers have applied OCT to minimize the expansion time with 
"frictionless conditions", i.e., taking an initial thermal equilibrium at one tem- 
perature into thermal equilibrium at another temperature in a cooling cycle, 
using real or imaginary bang-bang (piecewise constant or ramped) intermedi- 
ate trap frequencies, see e.g. Salamon et al. (2009); Hoffmann et al. (2011). 

3. 7 Other Applications 

Inverse engineering expansions using invariant theory or scaling laws have 
been applied in several contexts. For example, Choi, Onofrio and Sundaram 
(2011) discussed the possibility of achieving deep degeneracy of Fermi gases 
via sympathetic cooling by changing the trapping frequency of another species 
(the coolant) to keep constant the Lewis-Riesenfeld invariant. The identified 
advantages are the maximal heat capacity retained by the coolant due to the 
conservation of the number of atoms, and the preservation of its phase-space 
density in the nondegenerate regime where the specific heat retains its Dulong- 
Petit value. The limits of the approach are set by the transient excitation, 
that should be kept below some allowed threshold, and by the spreading of 
the cooling cloud which reduces the spatial overlap with the Fermionic cloud. 
The method is found to be quite robust with respect to broadband noise in 
the trapping frequency (Choi et al., 2012). 

Li et al. (2011) propose a scheme to cool down a mechanical resonator in a 
three-mirror cavity optomechanical system. The dynamics of the mechanical 
resonator and cavities is reduced to that of a time-dependent harmonic oscil- 
lator, whose effective frequency can be controlled through the optical driving 
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fields. A simpler harmonic system is studied in Zhang et al. (2012a), a charged 
mechanical resonator coupled to electrodes via Coulomb interaction controlled 
by bias gate voltages. 

Yuce (2012) designs, using scaling, fast frictionless expansions of an optical lat- 
tice with dynamically variable spacing (accordion lattice). Specifically, he con- 
siders the ID Hamiltonian H = p^/{2m) + V(t) cos (2A;z,x/A(t)) +ma;^(t)x^/2, 
where A is the scale parameter which goes from 1 at t = to c at t/ and 
the parabolic potential only acts during the expansion according to w^(t) = 
—hr^d'^Is./dt^. Decreasing the potential depth as V{t) = Vo/A^(t), and making 
the first and second derivatives of A vanish at the boundary times guarantee 
a frictionless expansion. In Ozcakmakli and Yuce (2012) the results are ex- 
tended to a continuously replenished BEC in a harmonic trap or in an optical 
lattice. 

Lau and James (2012) propose inverse engineering of the trap frequencies 
based on the Lewis- Riesenfeld invariants as part of the elementary operations 
necessary to implement a universal bosonic simulator using ions in separate 
traps. This method would allow to improve the accuracy and speed of conven- 
tional laser operations on ions which are limited by the Lamb-Dicke approxi- 
mation. 

Julia-Diaz et al. (2012) develop a method to produce highly coherent-spin- 
squeezed many-body states in bosonic Josephson junctions (BJJs). They start 
from the known mapping of the two-site Bose-Hubbard (BH) Hamiltonian 
to that of a single effective particle evolving according to a Schrodinger-like 
equation in Fock space. Since, for repulsive interactions, the effective potential 
in Fock space is nearly parabolic, the inversion protocols for shortcuts to 
adiabatic evolution in harmonic potentials may be applied to the many-body 
BH Hamiltonian. The procedure requires a good control of the time variation 
of the atom-atom scattering length during the desired period, a possibility 
now at hand in current experimental setups for internal BJJs. 



4 Transport 

The efficient transport of atoms and ions by moving the confining trap is 
a necessary fundamental requirement for many applications. These are for 
example quantum information processing in multiplexed trap arrays (Rowe 
et al., 2002; Reichle et al., 2006; Bowler et al., 2012) or quantum registers 
(Miroschnychenko et al., 2008); controlled translation from the production or 
cooling chamber to the interaction or manipulation zones; control of inter- 
action times and locations, e.g. in cavity QED experiments, quantum gates 
(Calarco et al., 2000) or metrology (Prestage et al., 1993); and velocity control 
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to stop (Schmidt et al., 2009; Chen et al., 2009) or launch atoms (Kuhr et al., 
2001). 

The transport should ideally be lossless, fast and "faithful" , i.e. the final state 
should be equal to the initial one apart from the translation and possibly phase 
factors. This is compatible with some transient excitation in the instantaneous 
basis at intermediate times. 

Many different experimental approaches have been implemented. Neutral atoms 
have been transported individually, as thermal atomic clouds, or condensates, 
using optical or magnetic traps. The magnetic traps can be displaced by mov- 
ing the coils mechanically, by time- varying currents in a lithographic pattern, 
or on a conveyor belt with permanent magnets (Lahaye et al., 2006). Op- 
tical traps can be used as optical tweezers whose focal point is translated 
by moving mechanically lenses (Convert et al., 2008a), and traveling lattices 
(conveyor belts) can be made with two counterpropagating beams slightly 
detuned. Mixed magneto-optical approaches are also possible. To transport 
ions, controlled time dependent voltages have been used in linear-trap based 
frequency standards (Prestage et al., 1993), and more recently in quantum 
information applications using multisegmented Paul traps (Huber et al., 2008; 
Walther et al., 2012; Bowler et al., 2012), or an array of Penning traps (Crick 
et al., 2010), also in 2D configurations (Blakestad et al., 2009). 

In general, a way to avoid spilling or excitation of the atoms is to perform 
a sufficiently slow (adiabatic) transport, but for many applications the total 
processing time is limited due to decoherence and an adiabatic transport may 
turn out to be too long. In the context of quantum information processing, 
transport could occupy most of the operation time of realistic algorithms, so 
"transport times" need to be minimized (Reichle et al., 2006; Huber et al., 
2008). There are in summary important reasons to reduce the transport time, 
and several theoretical and experimental works have studied ways to make fast 
transport also faithful (Convert et al., 2008a; Murphy et al., 2009; Masuda and 
Nakamura, 2010; Chen et al., 2010a; Torrontegui et al., 2011, 2012d). 



4-1 Invariant-based Shortcuts for Transport 

As done for expansions, shortcut techniques can be applied to perform fast 
atomic transport without final vibrational heating by combining dynamical 
invariants and inverse engineering. Two main scenarios can be handled in this 
way: shortcuts for the transport of a harmonic trap and shortcuts for the 
transport of an arbitrary trap. It is also possible to construct shortcuts for 
more complicated settings like atom stopping or launching, and combinations 
of transport and expansion of harmonic traps. 



29 



Transport of a rigid harmonic trap. - Suppose that a ID harmonic trap should 
be moved from qo{0) at time t = to d = qo{tf) at a time tf. The potential 
is V = Y^oi^ ~ Q'o(^))^ with fixed frequency. Comparing this to Eq. (7) this 
implies 

F = mu^qoit), uj{t) = Uq, U = 0. (54) 

Note that Eq. (8) plays no role here and Eq. (9) becomes the only relevant 
auxiliary equation, 

gc + Wo(9c-go) = 0, (55) 

where qc can be identified as a classical trajectory. This is the equation of a 
moving oscillator for which an analytical solution is known in both classical 
and quantum physics. From a classical mechanics point of view, the amplitude 
A of the oscillatory motion after transport is the modulus of the Fourier 
transform of the velocity profile associated with the trap trajectory (Convert 
et al., 2008a), 

^= |^[go](c^o)| (56) 
with J^lf] = /(t)e~*'^* dt. This Fourier formulation of the transport prob- 
lem allows for many enlightening analogies. For instance, A"^ is mathematically 
identical to the intensity profile for the far field Fraunhofer diffraction pattern 
of an object with a transmittance having the same shape as the velocity pro- 
file for the transport. An optimal transport condition is therefore equivalent 
to a dark fringe in the corresponding diffraction pattern. The optimization of 
the conditions under which a non adiabatic transport should be carried out 
with a rigid harmonic trap are thus equivalent to apodization problems in 
optics. If the velocity profile contains the repetition of a pattern one expects 
an interference-like effect, this would be, for instance, the case for a symmet- 
rical round trip transport as experimentally demonstrated in Convert et al. 
(2008a). 

Quantum mechanically, the wave function after transport reads 

= 4(, - ,M,t) exp ( ""'^ 7°">'* ) exp {l[ dt'm) . (57) 

where C = mg^/2 — mujQ^qc — qo{t)Y /2 is the Lagrangian associated with the 
equation of motion (55), and $ a wave function that coincides with the initial 
wave function at initial time and that evolves under the action of the static 
harmonic potential of angular frequency ojq located at g = ^0(0). Using the 
boundary conditions associated with the transport, one finds from Eq. (57) 
that an optimal transport for which the system starts in the ground state and 
ends up in the ground state of the displaced potential corresponds exactly to 
the classical criterion of a cancellation of the Fourier transform of the velocity 
profile, i.e. ^ = 0. 

Let us now address the application of invariant-based engineering. We first 
design an appropriate classical trajectory qcit) fulfilling the boundary condi- 
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tions g,(0) = go(0) = 0, = 0, qM = and q^tf) = qoitj) = d, q,{tf) = 

0) Qcitf) = 0, to ensure an evolution from the n-th state of the initial trap to 
the n-th state of the final trap. Then the trap motion trajectory qo(t) is de- 
duced via Eq. (55). Some variants are vertical transport with a gravity force, 
so that F = mu^qo — mg and Eq. (55) becomes qc + ^l^qc — qo) = —g, and 
stopping or launching processes (Torrontegui et al., 2011). 

A major concern in practice for all these applications is to keep the har- 
monic approximation valid. This may require an analysis of the actual po- 
tential and of the excitations taking place along the non-adiabatic trans- 
port process. Without such detailed analysis, the feasibility of the approach 
for a given transport objective set by the pair d,tf can be estimated by 
comparing lower excitation bounds (Torrontegui et al., 2012d). These are 
obtained using calculus of variations as we have discussed before for ex- 
pansions. Writing the expectation value of potential energy for a transport 
mode as {V(t)) = ^ [n + 1/2) + Ep^ the time average of Ep is bounded as 
'E^ > Qmd^/{tju^) (Torrontegui et al, 2011). 

This scaling should be compared to the milder dependence on tj"^ of the time- 
averaged transient energy in expansions (Chen and Muga, 2010). Chen et al. 
(2011b) have shown how to realise this bound by allowing the discontinuous 
acceleration of the trap at t = and t = tf and also finite jumps in the trap 
position. 

In Chen et al. (2011b), the invariant-based method is complemented by op- 
timal control theory. Since actual traps are not really harmonic, the relative 
displacement between the center of mass and the trap center is kept bounded 
as a constraint. The trajectories are then optimized according to different 
physical criteria: time minimization, (time-averaged) displacement minimiza- 
tion, and (time-averaged) transient energy minimization. The minimum time 
solution has a "bang-bang" form, and the minimum displacement solution is 
of "bang-off-bang" form. In this framework discontinuities in the acceleration 
(jc at the edge times and elsewhere are allowed. Physically this means that the 
trap may ideally jump suddenly over a finite distance, whereas the velocity qc 
and the trajectory qc remain always continuous. 

Transport of an arbitrary trap with compensating force. - In the second main 
scenario, the trap potential U{q — qo(t)) is arbitrary, and it is rigidly displaced 
along go(^)- Now, in Eq. (7), u = Uq = 0, F = rnqo, and qc in Eq. (9) may 
be identified with the transport function go- Inverse engineering in this case 
is based on designing the trap trajectory go (Torrontegui et al., 2011). In 
addition to U, there is a compensating linear potential term — mggo in H = 
p'^/2m — mq% + U{q — qo). The corresponding force compensates for the inertial 
force due to the trap motion in the rest frame of the trap, in such a way 
that the wave function in that frame is not modified up to a time dependent 
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global phase factor. This Hamiltonian was originally proposed by Masuda and 
Nakamura (2010) using the "fast-forward" scaling technique. Masuda (2012) 
has recently generalized this result for interacting, identical, spinless particles. 



4-2 Transport of a Bose-Einstein Condensate 



The two main scenarios of the previous subsection can be generalized for 
Bose-Einstein condensates (Torrontegui et al., 2012d). We first consider ID 
harmonic transport. For the GPE 



muj, 



2m dq^ ' 2 
the results of Sec. 4.1 motivate the ansatz 

{i i 
- {-fit + mq^q) - j^j^ dt' 

where x{(^) satisfies the stationary GPE 



+ ^{q-qo{t)y + 9imq.t)\' 



^(g,t), (58) 



m 



(59) 
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y2 ^ 1^:^^0 1,12 



a\' + U{a)+gM^)\' 



(60) 



The ansatz provides indeed a solution to Eq. (58) when qdt) satisfies Eq. (55). 
Inverse engineering gives the trap trajectory qo(t) from (55) after designing 
qdt), as for the linear dynamics. 

The inverse method can also be applied to anharmonic transport of conden- 
sates by means of a compensating force (Torrontegui et al., 2011). In either 
scenario this method does not require that tf satisfies any discretization condi- 
tion, as it occurs with other approaches (Torrontegui et al., 2012d), and tf can 
in principle be made as small as desired. In practice there are of course tech- 
nical and fundamental limitations (Torrontegui et al., 2011). Smaller values of 
tf increase the distance from the condensate to the trap center, and the effect 
of anharmonicity. There could be also geometrical constraints: for short tf, 
qo{t) could exceed the interval [0, d]. OCT combined with the inverse method, 
see below, provides a way to design trajectories taking these restrictions into 
account. 

Optimal control theory. — An OCT trajectory has been found when the cen- 
ter of the physical trap is kept inside a given range (e.g. inside the vacuum 
chamber), i.e. q^ < qo(t) < q^ (Torrontegui et al., 2012d). At the beginning 
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the trap is immediately set at the upper bound to accelerate the conden- 
sate as much as possible and at time ti the trap is moved to the lower bound 
to decelerate the condensate so as to leave it at rest aX tf. An important 
open question is to evaluate the effect of the approximate realization of the 
discontinuities found in the bang-bang solutions. 

Effect of Perturbations. -ToTTontegui et al. (2012d) also investigated the effect 
of anharmonicities when the harmonic transport protocol is applied. For a 
symmetrically perturbed potential V = u^m [{q — qoY + a{q — go)^] /2, the 
fidelity increases with increasing coupling constant gi, because of the increased 
width of the wavefunction. They also considered that the center of the physical 
trap is randomly perturbed with respect to qo{t). The fidelity at is found 
to be independent of d and the chosen qdt) and increases for shorter times tf 
and for smaller couplings gi, unlike the previous results. 



5 Internal State engineering 



Manipulating the internal state of a quantum system with time-dependent 
interacting fields is the basis of quantum information processing (Allen and 
Eberly, 1987; Vitanov et al., 2001; Bergmann et al., 1998) and many other 
fields. Two major routes are resonant pulses, and adiabatic methods such 
as "Rapid" Adiabatic Passage (RAP), Stimulated Raman Adiabatic Passage 
(STIRAP), and their variants. Simple fixed-area resonant pulses, such as a 
TT pulse, may be fast if intense enough, but they are also highly sensitive to 
variations in the pulse area, and to inhomogeneities in the sample (Allen and 
Eberly, 1987). Composite pulses provide an alternative to the single vr pulse, 
with some successful applications (Lewvitt, 1986; Collin et al., 2004; Torosov 
et al., 2011), but still they need an accurate control of pulse phase and in- 
tensity. In NMR, composite pulses are being superseded by adiabatic passage 
methods, which have also been very successful in laser cooling, chemical re- 
action dynamics, metrology, atom optics, interferometry, or cavity quantum 
electrodynamics. Adiabatic passage is robust versus parameter variations but 
slow. It is moreover prone to decoherence because of the effect of noise over the 
long times required. This motivates the search for fast and robust shortcuts, 
with respect to parameter variations and noise. 

Several methods to find STA have been put forward for two- and three-level 
atomic systems. Among them, methods that we have already discussed in Sec. 
2, like the transitionless driving, invariant-based engineering, or OCT. 
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5.1 Population Inversion in Two-level Systems 



Using the convention |1) = (^gj, |2) = y-^j, assume a two-level system with a 
Hamiltonian of the form 



2 \SjR(()+'in,(() A(i) 



In quantum optics it describes the semiclassical coupling of two atomic levels 
with a laser in a laser-adapted interaction picture, where Vtcif) = Qfi{t)+iQf(t) 
is the complex Rabi frequency and A(t) the time-dependent detuning between 
laser and transition frequencies. We will keep the language of the atom-laser 
interaction hereafter, but in other two-level systems, for example, in a spin- 
1/2 system or in a Bose- Einstein condensate on an accelerated optical lattice 
(Bason et al., 2012), Qc{t) and A(t) may correspond to different physical 
quantities. 

Initially at time t = 0, the atom is assumed to be in the ground state |1). The 
goal is to achieve a perfect population inversion such that at a time t = T 
the atom is in the excited state. For a it pulse the laser is on resonance, i.e. 
A(i) = for all t. If the Rabi frequency is chosen like Qc{t) = |fic(^)|e*", 
with a time-independent a, and such that Jq dt \^c{t)\ = vr, the population is 
inverted at time T. A simple example is the "flat" vr pulse with Qc{t) = e^'^ir/T. 

Adiabatic schemes provide another major route for population inversion. In 
the "Rapid Adiabatic Passage" technique the radiation is swept slowly through 
resonance. The term "rapid" here means that the frequency sweep is shorter 
than the life-time of spontaneous emission and other relaxation times. Many 
schemes corresponding to different functions A(t), and flR{t) are possible. The 
simplest is a Landau-Zener approach, with A linear in time, Qr constant and 
rii = 0. 

Transitionless shortcuts to adiabaticity.- If an adiabatic scheme is used and 
the adiabaticity condition ||fia| ^ 1^(^)1 (where f2 = y^A^ + f2/ = 0, 
and Qa = [^i?A — QrAI/Q"^) is not fulfilled, the inversion fails. We may still 
get an inversion (i.e. a shortcut) by applying a counterdiabatic field such that 
its maximum is not larger than the maximum of Qr (Chen et al., 2010c). 
The total Hamiltonian, see Sec. 2.2, for the transitionless shortcut protocol is 
(Demirplak and Rice, 2008; Berry, 2009; Chen et al., 2010c) 
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hi -A Qr- iQa 



(62) 



Invariant-based Shortcuts. -STA in two-level systems can be also found making 
use of Lewis- Riesenfeld invariants (Fasihi et al., 2012; Ruschhaupt et al., 2012). 
For Hq in Eq. (61), a dynamical invariant may be parameterized as 

h I cos(e(t)) sin(e(t))e— W\ 
lyt) = x/^ . h (63) 

^ \^sin(e(t))e*"W -cos(e(t)) y 

where /i is a constant with units of frequency to keep I{t) with dimensions of 
energy. From the invariance condition the functions 0(t) and a{t) must satisfy 

Q = Vtj cos a — VLr sin a, 

(64) 

a = —A{t) — cot (Qr cos a + f2/ sin a) . 
The eigenvectors of the invariant are 



\Mt)) 



cos(0/2)e-*"/2 
sin(0/2)e*"/2 



I 



V 



ia/2 



sin (0/2) e-' 



cos (0/2) e^°/2 



(65) 



with eigenvalues ±|yU. A general solution |\I/(t)) of the Schrodinger equation 
can be written as a linear combination 



m)) = c+e^^+(*)|0+(t)) +c_e^«-W|0_(t)), (66) 

where c± are complex, constant coefficients, and k± are the phases of Lewis 
and Riesenfeld (1969) introduced in Eq. (4). Let 7 = — 2k+ = 2k_, then 7 
must be a solution of 



7 = — ■ (cos a fiij + sin a fi/) . (67) 

sin 

Equivalently a solution of the Schrodinger equation |\E'(t)) may be designed 
with the same parameterization as above (|\E'(t))(\E'(t)| is a dynamical invari- 
ant.) and, by putting this ansatz into the Schrodinger equation, Eqs. (64) and 
(67) are found. 
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If flR(t), Qi{t) and A(t) are given, Eqs. (64) and (67) have to be solved to get 
6(t), a{t) and 7(t). A particular solution of the Schrodinger equation is then 
given by 



|^(t)) = |0+(t))e-^^W/^ (68) 

To find invariant-based shortcuts and inverse engineer the Hamiltonian 0(t), 
a{t), and 7(t) are fixed first, fulfilling the boundary conditions 6(0) = and 
0(T) = IT. The wave function (68) corresponds then to an atom in the ground 
state at t = and in the excited state at t = T, i.e. a perfect population 
inversion. Then, by inverting (64) and (67), 



Oii; = cosasin6 7 — sina 6, (69) 
f2/ = sin a sin 6 7 + cos a 6, (70) 
A = — COS67 — d. (71) 

There is much freedom in designing such a shortcut because the auxiliary 
functions Q{t), a{t) and 7(t) can be chosen arbitrarily except for the boundary 
conditions. 



5.2 Effect of Noise and Perturbations 



A key aspect to choose among the many possible shortcuts is their stability or 
robustness versus different perturbations. Ruschhaupt et al. (2012) have de- 
rived optimal invariant-based shortcut protocols, maximally stable concerning 
amplitude noise of the interaction and with respect to systematic errors. It 
turns out that the perturbations due to noise and systematic errors require 
different optimal protocols. 

Let the ideal, unperturbed Hamiltonian be the Ho{t) of Eq. (61). In Ruschhaupt 
et al. (2012), it is assumed that the errors affect Qf> and flj but not the de- 
tuning A, which, for an atom-laser realization of the two-level system is more 
easily controlled. 

For systematic errors, for example if different atoms at different positions are 
subjected to slightly different fields due to the Gaussian shape of the laser, the 
actual, experimentally implemented Hamiltonian is i^oi = Hq + f3Hi, where 
Hi{t) = HQ{t) and /3 is the amplitude of the systematic error. 

The second type of error considered in Ruschhaupt et al. (2012) is amplitude 
noise, which is assumed to affect Qr and flj independently with the same 
strength parameter A^. This is motivated by the assumption that two lasers 
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may be used to implement the two parts of the Rabi frequency. The final 
master equation describing systematic error and amplitude-noise error is 



7 ■ \ 2 



[H21, [H21, 



(72) 



where H2R{t) = Ho{t) 



A = nr = 



and H2i{t) = Ho{t) 



Before studying both types of error together it is fruitful to look at them 
separately. 

Amplitude-Noise Error. - If there is no systematic error (/3 = 0) and only an 
amplitude-noise error affecting the Rabi frequencies, a noise sensitivity can be 
defined as 



Qn ■= - 



1 d^P2 



2 9A2 



dPo 



A=0 



A=0 



where P2 is the probability to be in the excited state at final time T, i.e. 
P2 ~ 1 - QnX^. 

To find an invariant-based shortcut protocol maximally stable concerning am- 
plitude noise, it is first assumed that the unperturbed solution (68) satisfies 
6(0) = and 6(T) = vr. Using a perturbation approximation of the solution 
and keeping only terms up to (Ruschhaupt et al., 2012), 



1 , 
qN=-; at 
4 JO 



cos^ O + cos^ a sin^ 6) (m sin a — cos a6)^ 
+ (cos^ 6 + sin^ a sin^ 6)(m cos a + sinaO)^ 



(73) 



where m{t) = —7 sin 6. Minimizing the error sensitivity qn by Euler-Lagrange 
one gets that the optimal solutions satisfy (Ruschhaupt et al., 2012) a = nir/A, 
n odd, and 



(3 + cos(2e))e = sin(2e)(e)l 



(74) 



The corresponding Qr and Qi can be calculated from Eqs. (69) and (70). In 
this case, Qr = ±6/72 = ±ni and A(t) = 0. The optimal noise sensitivity 
value is = 1.82424/T < tt^/{4T) and the maximum of the Rabi frequency 
is QR(tf /2)tf ^ 2.70129. An approximate solution of Eq. (74) is given by 
Q{t) = 7it/T — ^ sin(27rt/T), with a noise sensitivity of = 1.82538/T. 

Systematic Error. — If there is no amplitude-noise error (A = 0) and only 
systematic error, a systematic error sensitivity is defined as 
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1 ^2^2 



dP2 

dm 



2 9/32 



where P2 is as before the probabihty to find the atom in the excited state at 
final time T. qs may be calculated with a perturbation approximation of the 
solution keeping only terms up to (3"^ (Ruschhaupt et al., 2012). 

To find an optimal scheme the invariant based technique is used again. The 
evolution of the unperturbed state can be parameterized as before, \ip{t)) (see 
Eq. (68)), with the boundary values G(0) = and Q{T) = 71. The expression 
for the systematic error sensitivity is now 



The optimal value is clearly qs = 0. An example of a class which fulfills 
g5 = is found by letting 7(t) = (20 — sin(26)) . If follows that qs = 
sin^ (nil) I (472^), so for = 1, 2, 3, ^5 = 0. There is still some freedom left, 
this allows further optimization concerning additional constraints. 

Systematic and amplitude-noise errors.- If both errors coexist the optimal 
schemes will depend on their relative importance. Ruschhaupt et al. (2012) 
examine numerically the behavior of different protocols. Fig. 1 shows that the 
different optimal schemes perform better than the other one depending on the 
dominance of one or the other type of error. 

Additional work is in required to extend the results in Ruschhaupt et al. (2012) 
to different types of noise and perturbations. Apart from the invariant-based 
approach, Lacour et al. (2008) have proposed robust trajectories in the adia- 
batic parameter space that maximize the population transfer for a two-level 
system subjected to dephasing. Also the robustness of the "parallel adiabatic 
passage" technique (keeping the eigenvalues of the Hamiltonian parallel (Dridi 
et al., 2009)) with respect to fluctuations of the phase, amplitude and pulse 
area was analyzed in Gulerin et al. (2011). 

5.3 Three-level Systems 

The transitionless driving for stimulated rapid adiabatic passage from level 1 
to level 3 in a lambda configuration with an intermediate state 2 making use 
of a pumping and a Stokes laser was studied in Unayan et al. (1997); Demir- 
plak and Rice (2003, 2005); Chen et al. (2010c). The fast-driving cd field 
connects levels |1) and |3). This implies in general a weak magnetic dipole 
transition, which limits the ability of the field to shorten the times. Invariant- 
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0.2 ' 0.00 



Fig. 1. (Color online) Probability P2 versus noise error and systematic 
error parameter; optimal systematic stability protocol (blue), optimal 
noise protocol (green). 

based engineering solves the problem by providing alternative shortcuts that 
do not couple directly levels |1) and |3) (Chen and Muga, 2012), as discussed 
below. It should be noted though that in an optical analogy of STA to engi- 
neer multimode waveguides all these schemes (with or without 1-3 coupling) 
may in principle be implemented (Lin et al., 2012; Tseng and Chen, 2012) by 
computer-generated holograms. In this analogy, based on the paraxial approx- 
imation, space plays the role of time so that the effect of the shortcuts is to 
shorten the length of the mode converters. 

In Chen and Muga (2012), using two lasers on resonance with the 1-2 and 
2-3 transitions, two single-mode protocols that make use of one eigenstate of 
the invariant are described. In these protocols full fidelity requires an infinite 
laser intensity, and shortening the time also implies an energy cost. The first 
protocol, based on simple sine and cosine functions for the pumping and Stokes 
lasers keeps the population of level 2 small. To achieve the same fidelity, less 
intensity is required in the second protocol, in which the intermediate level |2) 
is populated. The population of the intermediate level is usually problematic 
when its time decay scale is smaller than the process time. While this may 
be a serious drawback for an adiabatic slow process, it need not be for a fast 
shortcut. Protocols that populate level 2 may thus be considered as useful 
alternatives for certain systems and sufficiently short process times. 

In the previous two protocols the initial state is not exactly |1) to avoid a 
divergence in the Rabi frequency. A third multi-mode wave-function protocol 
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is also proposed in Chen and Muga (2012) using the same fields as for the 
first protocol but with an initial state which is simply the bare state |1). It 
provides a much less costly shortcut so exploring the multi-mode approach for 
this and other systems is an interesting task for future work. 

Inverse engineering of four-level systems has been considered in Giingordii et 
al. (2012), where a full Lie-algebraic classiffication and detailed construction 
of the dynamical invariants is provided. 



5.4 Spintronics 

Coherent spin manipulation in quantum dots is the key element in the state- 
of-the-art technology of spintronics. Ban et al. (2012) have considered the 
electric control of electron spin in a quantum dot formed in a two-dimensional 
electron gas confined by the material composition under a weak magnetic field, 
focusing on the spin fiip in the doublet of the lowest orbital state. The infiuence 
from higher orbital states can be taken into account by the Lowdin partition 
technique reducing the full Hamiltonian into an effective two-level one in which 
the matrix elements depend on electric field components. Using invariant- 
based inverse engineering the time-dependent electric fields are designed so to 
flip the spin rapidly and avoid decoherence effects. The results are stable with 
respect to environmental noise and the device-dependent noise and may open 
new possibilities for high-fidelity spin-based quantum information processing. 



5. 5 Experiments 

The counterdiabatic or transitionless approach described in Sees. 2.2, 2.4, and 
5.1 has been applied recently to invert the population of different two- level 
systems: 

In Bason et al. (2012) the effective two-level system is set as a condensate 
in the bands of an accelerated optical lattice Zenesini et al. (2009). Writting 
the Hamiltonian in Cartesian-like coordinates as H = Xax + YcTy + Zaz, X 
may be controlled by the trap depth, Z by the lattice acceleration Zenesini 
et al. (2009), and Y could in principle be implemented by a second shifted 
lattice. The counterdiabatic term in Eqs. (28) or (62) is of the form Yay whose 
realization is cumbersome in this setting. The alternative was to perform a 
unitary transformation that leads to the same final state modifying the original 
X and Z terms. This manipulation, discussed in Sec. 2.4. (see Eq. (31), was 
interpreted as a Z-rotation in Ibanez et al. (2012a), where it is compared 



40 



to the one based on the first-order superadiabatic cd-term H^^ . ^ Landau- 
Zener and a "tangent" protocol with a tangent function for Z unnafected by 
the rotation, are used as a reference, the later being found to be very robust 
versus a simulated variation of control parameters. 

In Zhang et al. (2012b) the two-level system is a single nitrogen vacancy cen- 
ter in diamond controlled by time dependent microwave fields. The reference 
process is a Landau-Zener transition, and the Yay cd-term is implemented 
by a field oscillating it /2 radians out of phase with respect to the field that 
provides the Xax term. As the maximal value of the total field amplitude is 
bounded, in this case to avoid undesired transitions, a "rapid-scan" approach 
is implemented to shorten the protocol time: the protocol is divided into a 
discrete set of time segments with varying phase and the time duration of 
each segment is adjusted so that the maximal amplitude allowed is applied. 



6 Wavepacket Splitting 

Splitting a wavefunction without excitating it is important in matter wave 
interferometry (Hohenester et al., 2007; Grond et al., 2009a,b; Pezze et al., 
2005). For linear waves, described by the Schrodinger equation, it is a peculiar 
operation, as adiabatic following is not robust but unstable with respect to a 
small external potential asymmetry (Gea-Banaloche, 2002; Torrontegui et al., 
2012b). The ground-state wavefunction "collapses" into the slightly lower well 
so that a very slow trap potential bifurcation fails to split the wave except 
for perfectly symmetrical potentials. A fast bifurcation with a rapidly growing 
separating potential succeeds to split the wave but at the price of a strong ex- 
citation. STA that speed up the adiabatic process along a non-adiabatic route 
overcome these problems (Torrontegui et al., 2012b). Numerical modelling 
shows that the wave splitting via shortcuts is significantly more stable than 
the adiabatic following with respect to asymmetric perturbations and avoids 
the final excitation. Specifically Torrontegui et al. (2012b) use the streamlined 
version (Torrontegui et al., 2012a) of the fast-forward technique of Masuda 
and Nakamura (2010), see Sec. 2.3, applied to Gross-Pitaevskii or Schrodinger 
equations after having found some obstacles to apply the invariants-based 
method (the eigenvectors of quadratic-in momentum invariants do not sat- 
isfy the required boundary conditions (Torrontegui et al., 2012a)), and the 
transitionless-driving algorithm (Demirplak and Rice, 2003) (because of dif- 
ficulties to implement in practice the counter-diabatic terms). The following 
discussion refers to the Schrodinger equation except for a final comment on 
the GPE. 

^ The use of the term "superadiabatic" in Bason et al. (2012) differs -is broader 
there- from the one in Sec. 2.4. 
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Fast-forward approach. — To apply the FF approach the density r{x,t) must 
be first designed. Assume the splitting of an initial single Gaussian /(x, 0) = 
g-x /(2aQ)^ where Oq is the width of the ground state for a harmonic oscillator 



with frequency u/h, oq = yh/{muj), into a final double Gaussian f{x,tf) = 
^-(x-xff/{2al) ^ g-(x+x^)V(2a2)_ rpj^g interpolation 

r(x,t) = z(t)[e-(---/W)'/(2-g) + e-(-+-/W)V(2ag)]^ ^75) 

where z{t) is a normalization function, generates simple F-shaped potentials. 
The conditions Xo(0) = io(^/) = are imposed, so r = at boundary times. 
In Torrontegui et al. (2012b) xq{s) = x/(3s^ — 2s^), where s = t/tf, is chosen 
for simplicity, and Eq. (25) is solved for the initial conditions to get the FF 
potential with Eq. (23). 

Effect of the perturbation. — The effects of an asymmetric perturbation may 
be studied with the potential Vx = Vpp + X9{x), where 9 is the step function 
and VpF the potential obtained via Eqs. (23), (25), and (75) with A = 0. The 
goal is to find a stable time-dependent porotocol that, even without knowing 
the value of A, is able to produce the split state. 

Moving two-mode model. — Static two-mode models have been used before to 
analyze splitting processes (Javanainen and Ivanov, 1999; Grond et al., 2009b; 
Aichmayr, 2010). Torrontegui et al. (2012b) consider instead a two- level model 
with moving left and right basis functions to provide analytical estimates and 
insight as a complement of the more detailed FF approach. 



Assume first the (symmetrical and orthogonal) moving left and right bare basis 



states \L{t)) = (°), \R{t)) = (l), and a corresponding two-mode Hamiltonian 




model 

m = l\ ^ ' ' \ , (76) 

where 6(t)/h is the tunneling rate (Javanainen and Ivanov, 1999; Grond et al., 
2009b) and A the energy difference between the two wells (Aichmayr, 2010). 
We may simply consider A constant through a given splitting process and 
equal to the perturbative parameter that defines the asymmetry. Thus, the 
instantaneous eigenvalues are E^(t) = ±|-y/A^ + S'^{t), and the normalized 
eigenstates 

|V'|(t)) = sin(f)|L(t))-cos(f)|i?(t)), 



\i;-{t)) = cos (f)|L(t))+sin(f)|i?(t)), 
where tana = S{t)/X defines the mixing angle. 



(77) 



When {\L{t)), \R{t))} are close enough initially (and 5(0) ^ A), the instan- 
taneous eigenstates of H are close to the symmetric ground state \ipo{0)) = 
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^(|L(0)) + |i?(0))) and the antisymmetric excited state |^(]"(0)) = ^(|L(0))- 
l-R(O))) of the single well. At two extreme regimes may be distinguished: 

i) For S{tf) > A the final eigenstates of H tend to IV^J (^/)) = -j^ilHtf)) ± 
\R{tf))) which correspond to the symmetric and antisymmetric splitting states. 

ii) For S{tf) <^ A the final eigenfunctions of H collapse and become right and 
left localized states: \ipx(tf)) = \L{tf)) and I'lpxitj)) = \R{tf)). Since S(tf) 
is set as a small number to avoid tunnelling in the final configuration, the 
transition from one to the other regime explains the collapse of the ground 
state function to one of the wells at small A ~ ^(^/)- 

Dynamics of the two-mode model. — In a moving-frame interaction-picture 
wave function = A'^ip^, where A = Z]/3=l,r l/5(^))(/5(0)l is the 

Schrodinger-picture wave function, ip^ obeys ihip"^ = {Ha — Ka)^^^, with 
Ha = A'^HA, and Ka = ihA^A. For real {x\R{t)) and {x\L{t)), the symmetry 
{x\R{t)) = {-x\L{t)) makes Ka = 0. 

We may invert Eq. (77) to write the bare states in terms of ground and ex- 
cited states, and get S(t) from Eq. (76). The actual dynamics is approximated 
by identifying \ijJo{t)) and -E'o'(^) with the instantaneous ground and excited 
states and energies of the unperturbed FF Hamiltonian. They are combined to 
compute the bare basis in coordinate representation, and with them the ma- 
trix elements {/3'\Hx\P). The dynamics in the moving frame for the two- mode 
Hamiltonian may then be solved. 

Sudden approximation. — The behaviour at low A may be understood with the 
sudden approximation (Messiah, 1999). Its validity requires (Messiah, 1999) 

tf < h/AlU, where ATU = \/(^(0)|:^'|^(0)) - {ij{0)\TUmoW and TU = 
Jq^ dt'HAit'). With \ip{0)) = \ipoiO)) the condition to apply the sudden 
approximation becomes A ^ In this regime the dynamical wave function 
ijj{tf) is not affected by the perturbation and becomes the ideal split state 
ipQitf), up to a phase factor. 

The previous results may be extended to weakly non-linear Bose- Einstein con- 
densates for which gi/a^X ^ 1. Otherwise the instability of adiabatic splitting 
with respect to perturbations is strongly suppressed by the compensating ef- 
fect of the non-linear term (Torrontegui et al., 2012b). Of course the shortcuts 
would still be useful if the time is to be reduced. 



7 Discussion 

We have presented an overview of recent work on shorcuts to adiabaticity 
(STA) covering a broad span of methods and physical systems. STA offer 
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many promissing research and application avenues with practical and fun- 
damental implications. Several pending tasks have been described along the 
text. We add here some more: To extend the set of basic physical limitations 
and laws for fast processes in specific operations, taking into account differ- 
ent constraints; To generate simple, viable shortcuts making systematic use of 
symmetries; To enhance robustness versus different types of noise and pertur- 
bations; To perform inverse-engineering with invariants beyond the quadratic- 
in-momentum family; To develop shortcuts for adiabatic computing, and in 
general for Hamiltonians that cannot be easily diagonalized, as in Nehrkorn 
et al. (2011); To design or supplement STA by optimal control theory meth- 
ods. We have seen some exaples but many other optimization problems await 
unexplored. 

Indeed STA open interesting prospects to improve or make realizable quantum 
information and technology operations, by implementing new fast and robust 
transport or expansion approaches, internal state manipulations, and cooling 
protocols; nuclear magnetic resonance is another field where developing ideal 
pulses may benefit from STA. STA could be also useful beyond single or many- 
body quantum systems, e.g. to build optical short-length mode converters or 
for designing mechanical operations with nanoparticles, mesoscopic, or macro- 
scopic objects. In classical mechanics there are many examples of adiabatic 
evolution that may be shortcut. The application of this concept to interacting 
classical gas manipulation remains also an open question. 

We have witnessed in a few years a surge of activity and applications that 
could have hardly been predicted. Researchers creativity will likely continue 
to surprise us in the stimulating crossroad of STA with new, unnexpected 
concepts and applications. 
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